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Abstract. The high redshift (z = 0.997) blazar 
B2 1308+326 was observed contemporaneously at x-ray, 
optical and radio wavelengths in June 1996. The x-ray ob- 
servations were performed with ASCA. The ASCA results 
were found to be consistent with reanalysed data from 
two earlier ROSAT observations. The combined ASCA 
and ROSAT data reveal an x-ray spectrum that is best fit 
by a broken power law with absorber model with photon 



= 1.63 



+0.10 



= 3.4_ 1 " 1 and Thard — j-"^_o.o9 
to* keV in the rest-frame of the 



spectral indices of r so f t = 
and a break energy at 1 . 1 
blazar. The break in the x-ray spectrum is interpreted, 
from the shape of the simultaneous broadband spectral 
energy distribution, to be the emerging importance of in- 
verse Compton (IC) emission which dominates the ASCA 
spectrum. The faint optical state reported for these obser- 
vations (my = 18.3 ± 0.25) is incompatible with the high 
synchrotron flux previously detected by ROSAT. The IC 
emission detected by both ROSAT and ASCA was not sig- 
nificantly affected by the large change in the synchrotron 
component. 

Mgn emission was detected with an equivalent width 
(Wa) of ~ 15A, significantly different from previously re- 
ported values. The small and variable Wa in B2 1308+326 
may be due to the highly variable continuum and not in- 
trinsically weak lines in the source. A lower limit on the 
Doppler boost factor calculated from the contemporane- 
ous data is consistent with expectations for highly po- 
larised quasars and higher than expected for BL Lacs. 
Absorption at a level of N H = 3.0±q;| x 10 20 cm~ 2 was 
detected which is in excess of the Galactic value of Nh = 
1.1 x 10 20 cm~ 2 , indicating the possible presence of a fore- 
ground absorber. A gravitational microlensing scenario 
cannot therefore be ruled out for this blazar. No signif- 



icant variability on timescales of hours was detected in 
the optical or x-ray data. 

B2 1308+326 could be a typical radio-selected BL Lac 
in terms of peak synchrotron frequency and optical and 
radio variability but its high bolometric luminosity, vari- 
able line emission and high Doppler boost factor make it 
appear more like a quasar than a BL Lac. It is suggested 
that B2 1308+326 be considered as the prototype of this 
class of composite source. 

Key words: BL-Lacertae objects: individual: 
B2 1308+326 quasars: individual: B2 1308+326 

X-rays: galaxies - Radio continuum: galaxies - gravita- 
tional lensing 



1. Introduction 

The classification of the hi gh redshift (z - 
B2 1308+326 is uncertain (iGabuzda et al. 



been designated a BL Lac (Stickel et al., 1991 



0.997 ) source 
1993|) . It has 
because 



of its almost featureless optical spec trum (Wills fc Wills 
1979 ), its strong optical variability ( Mufson et al., 1985 ) 
and the high degree of polarization of its optical contin- 
uum ( Angel fc Stockman, 1980[ ) . As one of the 1 Jy sample 
of BL Lacs, B2 1308+326 is classified as a Radio-Selected 



BL Lac (RBL) ( [Stickel et al., 1993| ). 

However, B2 1308+326 also exhibits quasar-like prop- 
erties. VLBI polarisation images show that the polarised 
flux from the inner part of the jet is nearly perpendic- 
ular to the direction of the jet which is a characteristic 



property of quasars not BL Lacs (Gabuzda et al., 1993) 
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The degree of core polarisation from VLBI observations of 
B2 1308+326 is somewhat higher than usual for quasars 
|and at the lower end of the range for BL Lacs jGabuzda 
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et al., 1993). Measurements of the apparent speed of 
knots emerging from the core are superluminal (3.6/i -1 c, 
8h~ 1 c, 21h~ 1 c) and above the typical speeds observed in 
BL Lacs (pabuzdact al., 1993|). B2 1308+326 also shows 



an optical- UV excess above the IR extrapolation (Brown 
et al., 1989), and its radio power and morphology are 



FRII (Fanaroff fc Riley, 1974, type II) in nature (Koll 



gaard 3t al., 1992). The FRII classification is not that 



unusual among BL Lacs (Kollgaard et al., 1996), but an 



FR I morp hology would strongly militate against its being 
a quasar flUrry et al., 2000|; JKoTlgaard et al., 1990|) . The 
large bolo metric (Sambruna et al., 1996 ) and 5 GHz lu- 
minosities ( Kollgaard et al., 1992 ) are further indications 
that it could be regarded as a quasar. The term "blazar" 
is particularly useful for sources like this, encompassing as 
it does flat spectrum radio quasars (FSRQs) and BL Lacs. 
It must be noted that the terms highly polarised quasars 
(HPQs), optically violently variable (OVV) quasars, FS- 
RQs (and more recently, quasar-hosted blazar) are used to 
refer to blazar subclasses. The significant distinction here 
is between quasar-type and BL Lac-type blazars and the 
term FSRQ is used except in cases where another classifi- 
cation occurs in a work cited in the text. 

Two main approaches have proven to be effective in 
constraining models of BL Lacs; monitoring their variabil- 
ity (e.g. Georganopoulos & Marscher 1998a), and studying 
their broadband spectral energy distributions (SEDs) (e.g. 
Sambruna et al. 1996, Kubo et al. 1998, Sambruna et al. 
1999). Rapid flux variations and changes in the shape of 
the SEDs of blazars are common (e.g. Pian et al. 1998, Ma- 
comb et al. 1995, Chiaberge & Ghisellini 1999). To place 
adequate constraints on the synchrotron emission (e.g. the 
Doppler boost factor of the putative jet and the angle be- 
tween such a jet and the line of sight) in these sources 
requires simultaneous observations over a wide range of 
frequencies. 

It has been suggested that some sources with BL Lac 
characteristics are actually gravitationally microlensed 
quasars flNottale, 19"86| ; pstriker fc Vietri, 1990|). The bes t 
candidates include AO 0235+164 (|Rabbette et al., 19961 ), 



PKS 0537-441 ([Romero et al., 1999| ) and MS 0205.7+3509 
( |Watson et al., 1999| ). B2 1308+326 has characteristics in- 
termediate between BL Lacs and quasars. Gabuzda et al. 
( |1993| ) have therefore suggested that B2 1308+326 may be 
gravitationally microlensed. Observation of a foreground 
galaxy would lend support to the argument that it is in- 
deed gravitationally lensed. However high resolution imag- 
ing with the HST WFPC2 ( |Urry et al., 1999[ ) did not 
confirm the detection of spatially extended emission re- 
ported by Stickel et al. (1993). The observations were 
consistent with a point source with an upper limit of 
mi > 26magarcsec -2 for any surrounding galaxy. 

In the superluminal micro-lensing scenario, the time 
scale for rapid variability can constrain source parame- 
ters (e.g. Lorentz factor and source diameter, McBreen 
& Metcalfe 1987, Gopal-Krishna & Subramanian 1991). 



Such microlensing variations should be achromatic and 
therefore detectable with simultaneous multi-wavelength 
observations. 

Simultaneous x-ray, optical and radio observations 
were carried out in order to constrain the nature of the 
emission in B2 1308+326. These observations and the data 
reduction procedures are described in Sect. 2. Results are 
presented in Sects. 3 and 4 and are discussed in Sect. 5 
with conclusions in Sect. 6. Hq = 75kms~ 1 Mpc~ 1 and 
<7o = 0.5 are assumed throughout. 

2. Observations and data reduction 

2.1. X-ray observations 

Observations of B2 1308+326 were made with the ASCA 



satellite (Tanaka et al., 1994) between 15:26 and 23:41 UT 



on 10 June 1996, and between 23:37 UT on 11 June 1996 
and 09:10 UT on 12 June 1996. There are four instru- 
ments on board ASCA, two Gas Imaging Spectrometers 
(GIS-S2 and GIS-S3) and two Solid-state Imaging Spec- 
trometers (SIS-S0 and SIS-SI). The instruments have a 
well-calibrated energy range of 0.7-10.0 keV for the GIS 
and 0.5-10.0 keV for the SIS and only data in these en- 
ergy ranges were used for spectral fitting. The total ex- 
posure time per instrument was approximately 40 ks. The 
SIS data were recorded in FAINT data mode and con- 



verted to BRIGHT-2 mode (Inoue, 1993) having corrected 
for the Residual Dark Distribution (RDD) effect with the 
FTOOLS script 'Correctrdd'. 

Data from the SIS and GIS instruments were screened 
in a standard way, and were then reduced using the 
FTOOLS applications. An extraction region was defined 
around the source, with an aperture of 4' for the SIS de- 
tectors and 6' for the GIS. This spatial region was used 
to extract the source plus background counts. An esti- 
mate of the background was derived by taking an annulus 
around the extraction aperture and extracting counts from 
this region. An alternative background estimate was de- 
rived from separate dark-sky observations made by ASCA 
and results obtained using both types of background es- 
timates were found to be consistent. The annulus back- 
ground was subtracted from the extracted source plus 
background counts to obtain a source spectrum. The re- 
sulting spectra were rebinned to have a minimum of 20 
counts per bin in order to ensure the validity of the Gaus- 
sian approximation. 

Two observations of the source were made with the 



ROSA T PSPC dPfcffcrmann fc Briel, 1986|) on 2 3 June 

1991 dComastri ct al., 1995| ; |Lamer et al., 1996J ) and 3 
June 1992. These data were reduced in a standard way 
using the FTOOLS software. The response matrix from 
ROSAT AO-1 was used to analyse the data from the 23 
June 1991 observation, and that from AO-2 for the 3 June 

1992 observation. Only the well-calibrated data between 
0.1-2.4keV were used. 
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2.2. Optical and radio observations 



Optica J Johnson V-band observations were carried out un- 



(Table |l|) and to model B, a broken power-law with ab- 
sorber (Table ph. The /-test (Bevington & Robinson 



der ph otometric conditions with the 1.23 m telescope at 
Calar Alto Observatory, Spain on the nights of 10 and 11 
June 1996 with typical exposure times of 900 s. Standard 
CCD aperture photometry techniques were used to reduce 
the data. Absolute calibration of the source was derived 
using V-band differential photometry on a calibrated ref- 
erence star within the CCD field of view (star C from 
Fiorucci et al. 1998). Optical V-band observations made 
two weeks later with the 1 m JKT telescope at La Palma 
Observatory showed the source to be of the same mag- 
nitude, within the error of 0.25 mag. The relative uncer- 
tainty (0.08 mag.) in the optical magnitude was brought 
below the absolute calibration uncertainty (0.25 mag.), by 
using differential photometry techniques ( Ftabbette et al 
199SD 



1992|) was applied to models A and B to establish whether 



Observations were carried out under non-photometric 
conditions at the 2.12 m telescope of the Observatorio 
"Guillcrmo Haro" at Cananea, Son., Mexico on June 16, 
1996 using the Landessternwarte Faint Object Spectro- 
graph and Camera (LFOSC) in ^5.5 A/pixel mode. The 
standard IRAF utility package distributed by the National 
Optical Astronomy Observatories was used to reduce the 
data. The spectral images were processed to obtain one- 
dimensional spectra of the object and the sky by co-adding 
the counts on each channel along the slit, through user de- 
fined apertures. 

The wavelength calibration was performed using the 
identified lines of the comparison arcs in order to obtain a 
dispersion solution. With the derived coefficients the wave- 
length calibration of the spectra was then carried out. Due 
to the non-photometric conditions during these spectral 
observations, the spectrophotometric calibration proved 
to be unreliable. 

The results of radio observations at 22 GHz and 
37 GHz performed at the Metsahovi Radio Observatory, 
Finland for the period 1980-2000 are included here. The 
observing and data reduction procedures are described in 
Terasranta et al. Q1998J ). 



3. X-ray Spectral Fits 

The x-ray data from both ASCA observations were com- 
pared for variability above a significance level of 2 a. No 
variability was detected and the data were therefore coad- 
ded per detector for further analysis. 

Data from the four ASCA instruments were fit in- 
dividually using the Levenberg-Marquardt algorithm in 
XSPEC. Model A, a power-law with an absorber as de- 
fined by Morrison & McCammon ( 1983| ) was fit to these 
data and yielded acceptable fits and consistent results (Ta- 
We0). 

Since the ASCA instruments showed no significant 
spectral differences, the binned spectra from all four detec- 
tors were fit simultaneously. The data were fit to model A 



the fit was significantly improved by the addition of the 
two extra terms in model B (Table ph. Model B did not 
yield a significant improvement over model A in the fit to 
the ASCA data alone (Table |). 

The ROSAT da ta fro m 23 June 1991 were published 
by Comastri et al. ( 1995| ). The results of a single power- 
law with absorber fit (model A) to this ROSAT data were 
found to be compatible with and had a similar reduced \ 2 
statistic, xt (Table [j]) to that obtained by Comastri et al. 



(1995) for the same data and model. Although model A 
provides an acceptable fit (xt = 1-13) to the ROSAT data 
of 23 June 1991 (Rl), the /-test indicates that a broken 
power-law is a significantly better fit (Table |||). The spec- 
tral break is clearly detected in this data alone (Fig. m. 
Data from the ROSAT observation of 3 June 1992 (R2) is 
acceptably fit by model A; there is no significant improve- 
ment in the fit to this data using model B (Table y) . 

The energy of the spectral break is at the low energy 
end of the SIS detector range and outside that of the CIS. 
The ROSAT observation of 3 June 1992 unfortunately 
does not constrain the models very well. It is not sur- 
prising therefore that there is no significant improvement 
in the fit for either of these observations using a broken 
power-law model. Data from both ROSAT observations 
were fit with model B and the results were found to be 
consistent with the results from the fit of model B to the 
ASCA data (Table ||). No variation was found at the level 
of 95% significance, in the 0.5-2.4keV flux between any 
of the observations. 

The consistent results indicate that the x-ray spectrum 
had not altered substantially between the ROSAT and 
ASCA observations (see the next section for further dis- 
cussion on this point). The combined ROSAT and ASCA 
data were therefore fit simultaneously, allowing the nor- 
malization parameter to vary independently between the 
two ROSAT observations and the combined ASCA data, 
but constraining the other parameters to be the single 
best-fit for all datasets. The flux normalizations were left 
free to vary independently to allow for any residual sys- 
tematic differences between the instruments. 

Model B was the best fitting model to the combined 
ASCA and ROSAT data and was used for all further 
analysis. Although a discrepancy of Ar ~ 0.4 has been 
noted in the spectral indices reported for observations with 
ASCA SIS and ROSAT PSPC Jlwasawa et al., 1999| ), the 
detection of the spectral break does not depend on dif- 
ferences between the ROSAT and ASCA spectral indices 
since the break is clearly detected in the Rl data alone. 
The results of a comparison between models A and B for 
all the x-ray observations are summarised in Table 0. 
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Table 1. Comparison of fits (using model A, a simple 
power-law plus absorber) to data from: each ASCA detec- 
tor individually (SO, SI, S2, S3), from the combined ASCA 
detectors (A), from the ROSAT observations of 23 June 
1991 (Rl) and 3 June 1992 (R2) and from both ASCA 
and ROSAT, fit simultaneously (A+R) . Nh is the equiv- 
alent hydrogen column, T is the photon index, Fik v is 
the flux in photons/keV/cm 2 /s at IkeV, x 2 /DOF is the 
value of x 2 for that fit with DOF degrees of freedom and 
N. H. Prob. is the probability of the x 2 obtained being at 
least as large as it is if the data come from a parent dis- 
tribution described by the model. Ranges in parentheses 
are 90% confidence intervals for 1 parameter of interest. 



(10 



N H 

20 „„ 



'I 



(xlO" 



X7DOF 
N. H. Prob. 



SO 



SI 



S2 



S3 



Rl 



R2 



A+R 



0.0 
(0.0-11.6) 

0.4 
(0.0-14.6) 

24.5 
(0.0-102.4) 

0.0 
(0.0-22.3) 

0.0 
(0.0-2.5) 

1.35 
(0.84-1.94) 

0.81 
(0.00-2.15) 

1.36 

(1.14-1.58) 



1.45 
(1.28-1.71) 

1.38 
(1.18-1.74) 

2.21 
(1.59-3.24) 

1.83 
(1.52-2.23) 

1.75 
(1.66-1.84) 

2.01 
(1.76-2.27) 

1.75 
(1.29-2.30) 

1.88 
(1.80-1.96) 



1.94 
(1.67-2.55) 

1.82 
(1.52-2.60) 

2.94 
(0-9.20) 

2.42 
(1.82-3.59) 

1.99 
(1.87-2.16) 

2.19 
(1.99-2.39) 

1.89 
(1.61-2.17) 

2.15 
(2.05-2.25) 



46/45 
0.45 

29/38 
0.84 

37/49 
0.90 

67/61 
0.27 

205/202 
0.43 

34/30 
0.29 

10/12 
0.66 

250/247 
0.43 



4. Results 

The deconvolved data from all instruments with the best- 
fit model and the data/model ratio are plotted in Fig. [TJ. 
Confidence contours for the best-fit model for the two 
parameters of interest, equivalent Hydrogen absorbing col- 
umn (Nh) and the spectral break energy, are shown in 
Fig. 0. The soft x-ray absorption at 3.0 x 10 20 cm -2 is 
above the Galactic level of 1.1 x 10 20 cm -2 given by Dickey 



& Lockman ( |1990| ). 

The rest-frame SED is plotted in Fig. 0. The break 
in the ROSAT spectrum is clearly evident at 1.1 keV 
(2.8 x 10 17 Hz). It is generally accepted that the two broad 
emission components that characterise blazar SEDs are 
due to synchrotron emission at lower energies and in- 



verse Compton (IC) emission at higher energies (Costa- 
mante et al., 2000]; jFossati et al., 2000 ; Sambruna et al. 
1996 |; |yian et al., 1996Q The SED of B2 1308+326 shows 
the synchrotron component peaking between 10 12 Hz and 
10 14 Hz with the inverse Compton emission dominating 



Table 2. Results from best- fit broken power-law plus ab- 
sorber model (B), to data from different observations; Rl 
and R2 are the ROSAT observations of 23 June 1991 and 
3 June 1992 respectively; A represents the simultaneous 
fit to the data from the ASCA observations of 10-12 June 
1996; A+R is a simultaneous fit to all the ROSAT and 
ASCA data. Ranges in parentheses are 90% confidence 
intervals for 1 parameter of interest. Tl and T2 are the 
photon spectral indices below and above the break energy 
respectively. 



N H 
(10 20 cm- 2 ) 



n 

(soft) 



Break 
(keV) 



T2 

(hard) 



Rl 



"A 



A+R 



8.4 
(6.35-8.69) 

7.1 
(5.1-7.6) 

2.78 
(0.0-9.7) 

3.0 



> 10.0 

(3.04-) 
10.0 

(") 
10.0 

(-) 

3.4 



0.38 
(0.35-0.43) 

0.33 
(0.28-0.42) 

0.45 
(0.21-0.56) 

0.57 



2.1 
(1.7-2.4) 

2.32 
(1.58-2.49) 

1.62 
(1.54-1.77) 

1.63 



(2.4-5.3) (2.3-8.5) (0.36-0.77) (1.54-1.73) 



a n is unconstrained for the ASCA and R2 data, since 
model B is not a significantly better fit than model A for these 
data alone (see Table [3J) . 

Table 3. Fit comparisons for models A, power-law 
plus absorber and B, broken power-law plus absorber. 
N. H. Prob. is the probability of the % 2 obtained being at 
least as large as it is if the data come from a parent distri- 
bution described by the model. The /-test Prob. refers to 
the probability derived from the /-test of the hypothesis 
that an additional component (i.e. spectral break) yields 
an improved description of the data. Rl and R2 are the 
ROSAT observations of 23 June 1991 and 3 June 1992 re- 
spectively; A represents the ASCA observations of 10-12 
June 1996; A+R is a simultaneous fit to all the ROSAT 
and ASCA data. The spectral break is not significant in 
the ASCA or R2 data individually, but is significant in Rl 
alone and in the A+R simultaneous fit. 



Obs. 


A: absxpow 
N.H. Prob. ^ 


B: absxbknpow 
N.H. Prob. J& 


/-test 
Prob. 


Rl 
R2 

A 
A+R 


0.286 34/30 
0.657 10/12 
0.431 205/202 
0.430 250/247 


0.707 24/28 
0.657 7.7/10 
0.433 203/200 
0.621 238/245 


0.994 
0.651 
0.645 
0.998 



by ~ 2.8 x 10 17 Hz. The parabolic fit to the contem- 
poraneous data indicates that the peak synchrotron fre- 
quency is ~ 10 12 Hz (Fig. 0). Although no variation in 
the 0.5-2.4 keV spectra was detected between the observa- 
tions made with ROSAT and ASCA, examining the SED 
(Fig. pf) leads to the conclusion that the ROSAT syn- 
chrotron flux is incompatible with the radio and optical 
data obtained contemporaneously with ASCA. The effect 
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Fig. 1. Best fit (observer frame) ExFe spectrum to joint ASCA and ROSAT data with data/ model ratio. The break 
at 0.57kcV is indicated. Data from the ASCA instruments are plotted in grey (SO 'V', SI '▼', S2 'A' and S3 'A'). 
ROSAT observations are plotted in black (Rl (23 June 1991) '•', R2 (3 June 1992) 'o'.) 



of the absence of the synchrotron flux on the Rl spectrum comb ct al., 1995 ). Their model allows a peak synchrotron 



was quantified by extrapolating the best-fit ROSAT syn- 
chrotron power-law past the spectral break to 2.4 keV. The 
0.5-2.4 keV flux of this synchrotron power-law alone was 
subtracted from the composite synchrotron-IC Rl 0.5- 
2.4 keV flux. The resulting (IC) flux was still fully con- 
sistent with the ASCA 0.5-2. 4 keV flux. This indicates 
that the synchrotron-IC break occurred at lower energies 
during the ASCA observation than during Rl and that 
the ASCA spectrum is almost entirely dominated by an 
IC component that has not varied substantially at these 
energies despite the large change in the synchrotron flux. 
Chiaberge & Ghisellini ( |1999[ ) have modelled synchrotron- 
self-Compton emission in blazars and successfully repro- 
duced the SED of Mkn421 in different flux states (|Ma- 



variation of nearly two orders of magnitude that leaves 
the low-frequency end of the IC component largely un- 
affected, consistent with the SED observed here and in 



Mkn421 flMacomb et al., 1995[ ). Major variability of the 
IC component would however be expected at 7-ray ener- 
gies because of the large change in the synchrotron emis- 
sion. Observations at 7-ray energies can strongly con- 
strain synchrotron-IC models (Maraschi ct al., 1999). A 



parabolic fit to the highest archival fluxes in each wave- 
band and the ROSAT data below the spectral break yields 
a peak synchrotron frequency of ~ 10 14 Hz. This indicates 
that the peak synchrotron frequency of B2 1308+326 may 
decrease with flux, an important indicator in deciding the 
nature of the inverse Compton emission in this source 



D. Watson et al.: ASCA and other contemporaneous observations of the blazar B2 1308+326 




0.04 0.06 0.08 

N H (x 10 22 cm ) 

Fig. 2. 68%, 90% and 99% confidence contours for two 
parameters of interest; equivalent hydrogenic absorbing 
column (Nh) and spectral break energy, both in the ob- 
server frame. The dashed line represents the equivalent 
Hydrogen column density level of the Galaxy in the direc- 
tion of B2 1308+326 ([Dickey fc Lockman, 1990|). 



(Bottcher, 1999). This type of spectral break and the pos- 
sibility of increasing synchrotron peak frequency with in- 
creasing flux have also been observed with BeppoSAX in 
the blazar S5 0716+714 (|Giommi et al. 



1999) 



The ASCA lightcurve yields no evidence for variabil- 
ity above 7.5 x 10~ 14 ergcm~ 2 s _1 on timescales of hours. 
V-band optical photometry from both Calar Alto observa- 
tions yields a magnitude of 18.3 + 0.25. A limit on the vari- 
ability of the source of 0.08 mag. was found using differen- 
tial photometry. Therefore, no constraints can be placed 
on the supcrluminal gravitational lensing scenario in this 
AGN because of the absence of variability. 

The optical spectrum (Fig. ||) observed on 16 June 
1996 is intermediate between the extreme spectra previ- 
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Fig. 3. Rest-frame, v¥ v spectral energy distribution of 
B2 1308+326. The x-ray data has been corrected for ab- 
sorption using the best-fit parameters (Table ||). The 
break at 2.8 x 10 17 Hz is clearly seen in the middle of 
the x-ray data. Data contemporaneous with the ASCA 
observation are plotted in black. Non-simultaneous data 
are plotted in grey. Archive data are from NEDq Two 
parabolic fits to the synchrotron component are plotted, 
the black one is the best fit parabola to the contemporane- 
ous data below the break, the grey parabola is a fit to the 
highest archival flux values in each band (from NED) and 
the ROSAT data below the break. The shift in frequency 
of the synchrotron peak is apparent. 



first discovered, recent work has focused on determining 
the physical distinctions between the t wo types ( Padovani| 



The eq uivalent width of the Mg II line was ~15 A. The ab 



solute flux calibration for this spectrum is unreliable due 
to non-photometric conditions during the observation; it 
was therefore not possible to obtain a continuum level or 
line fluxes from the ASCA observations. 

The results of the radio observations (Fig. 5) over 20 
years indicate that B2 1308+326 was in an intermediate 
radio flux state during these observations. 

5. Discussion 

Although BL Lacs were originally classified as XBLs or 
RBLs depending on the energy band in which they were 

1 The NASA/IPAC Extragalactic Database (NED) is oper- 
ated by the Jet Propulsion Laboratory, California Institute of 
Technology, under contract with the National Aeronautics and 
Space Administration. 



[Fossati et al., 1997 ; Gcorganopoulos & 
Marscher, 1998b| ) and their relation to the other subclass 
of blazar, FSRQs (Sambruna et al., 2000J; |Sambruna, 1997; 



Fossati et al., 199£; phisellini et al., 1998 ). The spectral 
energy distributions of XBLs and RBLs differ markedly, 
with the peak in synchrotron output generally occurring 
in the x-ray range for the former and the infrared range for 
the latter flPadovani fc Giommi, 1995| ). This property has 
allowed BL Lacs to be classified as low-frequency peaked 
BL Lacs (LBLs, roughly the same sample as RBLs) and 
high-frequency peaked BL Lacs (HBLs, approximate to 
XBLs). Such differences may arise due to changes in the 
dominant radiation mechanism occurring as a function of 
both the angle between the line of sight and the bulk ve- 
locity (phisellini fc Maraschi, 1989), although it has been 



demonstrated, assuming certain conditions, that the com- 
plete range of properties observed in blazars are difficult 
to account for solely by a change in orientation of the jet 
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Fig. 4. Rcdshift (z — 0.997) corrected optical spectrum. 
The absolute flux calibration is unreliable due to non- 
photometric conditions during the observation, but is in- 
cluded for completeness. The Mgll emission feature is 
identified. 



flBambruna ct al., 1996 ). On the other hand, XBLs and 
RBLs have different variability timescales, a fact that is 
very naturally explained with the orientation hypothesis 
( Hcidt fc Wagner, 1998 ). The recent RGB sample revealed 
BL Lacs intermediate between XBLs and RBLs, indicat- 
ing that there is no bi-modal distribution in BL Lacs and 
that they span peak synchrotron frequencies continuously 
from IR to x-ray wavelengths ( Laurcnt-Muchlciscn ct al 
19990 

One hypothesis connecting FSRQs and BL Lacs is 
that t hey arc high and low luminosity radio galaxies re - 



spectively that are highly- relativistically beamed towards 
the observer (Padovani, 1992). A more recent theory con- 
nects FSRQs, LBLs and HBLs in a sequence of increasing 
peak s ynchrotron frequencies and decreasing bolomctric 



luminosity ( Fossati et al., 1998 ; Bambruna et al., 2000| ); 
however, roughly 2 5% of FSRQs in the D eep X-ray Ra- 
dio Blazar Survey ( Perlman et al., 1998 ) have HBL-like 



radio-x-ray broadband spectral indices (artx), indicating 
that the synchrotron peak frequencies (t'peak) for FSRQs 
reach energies at least as high as those of HBLs, since 
c*rx is strongly correlated with i/ pca k ( Fossati et al., 1998 , 
and references therein). In this theory, blazars with strong 
thermal emission or strong emission lines (FSRQs) should 
have soft x-ray spectra dominated by inverse Compton 
(IC) emission produced by jet electrons scattering pho- 
tons produced externally to the jet (external Compton, 
EC model). Blazars with intrinsically weak emission lines 
(BL Lacs) are expected to have their IC component dom- 
inated by jet electrons scattering photons produced in the 
jet via synchrotron emission (synchrotron self-Compton, 
SSC model). Bottcher ( |1999| ) has postulated that EC- 
dominated spectra have peak synchrotron frequencies that 
decrease with increases in flux. The peak synchrotron fre- 
quency of B2 1308+326 appears to decrease with decreas- 
ing flux. While this synchrotron peak shift behaviour is 
contra: y to the trend found in samples of blazars ( Kubo 
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Fig. 5. 1980-2000 radio lightcurve at 22 GHz (top) and 
37 GHz (bottom) obtained at the Metsahovi Radio Obser- 
vatory. The times of the ASCA and ROSAT (Rl and R2) 
observa ;ions are indicated. 



et al., 1998; Fossati ct al., 1998 ) it is not at all unusual in 



flares detected in individual sources (e.g. Catanese et al. 
1997). Without IR photometry it is not possible to deter- 
mine precisely the peak synchrotron frequency at the time 
of these observations. The optical observations were ob- 
tained simultaneous with, and two weeks after, the ASCA 
observations and yielded compatible results on both oc- 
casions. The radio data nearest the ASCA observations 
were taken on 4 June 1996 at 37 GHz, and on 5 and 20 
June 1996 at 22 GHz. The radio variability during the sur- 
rounding months was less than 30%. The spectrum would 
have to deviate significantly from a parabolic shape in or- 
der to peak above 10 14 Hz and still manage to drop as 
low as 18.3 mag. in the optical. Such a spectral shape can- 
not be ruled out without accurate modelling of the syn- 
chrotron component, however it is very likely that the syn- 
chrotron peak did occur below 10 14 Hz in the faint state 
observed here, indicating that the synchrotron peak fre- 
quency decreases with decreasing flux and that SSC is the 



more likely emission mechanism in this blazar ( Bottcher 



1999 ). Simultaneous GeV observations are required to de- 
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termine if there is excess IC emission resulting from the 
EC component. Such observations would be useful for test- 
ing the model proposed for PKS 0548+134 (Mukherjee 
et al., }999D and the trend detected by Kubo et al. (|1998|) 



in a sample of 18 blazars. 

The dominance of the IC component in B2 1308+326 
above 2.8 x 10 17 Hz is evident from Fig. || The relatively 
low energy of this break is consistent with expectations 
of LBLs ( Giommi ct al., 1995J ), but somewhat exceeds ex- 
pectations for FSRQs in the blazar paradigm described by 



Sambruna et al. (2000). The contemporaneous nature of 



the observations rules out distortion of the spectrum due 
to source variability. 

The primary reason for classifying B2 1308+326 as a 
BL Lac rather than a FSRQ is its largely featureless opti- 
cal continuum which is one of the defining features of the 



BL La c class, with rest- frame Wa < 5 A for any line (Sam- 
bruna ^t al., 199q ) . This criterion has been j ustifiably crit- 



icise d as arbitrary (|Kollgaard ct al., 1992 ; Marcha et al 



1996; Scarpa fc Falomo, 1997 ) partly on the basis that 
some sources (such as B2 1308+326) have Wa that vary 
substantially across this boundary. Nevertheless, weak or 
absent emission lines is a distinguishing property of these 
sources (Padovani, 1992). The rest-frame Wa of the Mgn 
line in B2 1308+326 detected by Miller et al. ( |1978|) wa s 
< 5 A and no lines were observed by Wills & Wills ( 1979Q . 
However, an Mgn Wa = 18.6 A in the rest-frame was ob- 
served in March 1989 by Stickel et al. (1993) in a very 
low luminosity state m.R = 19.0 mag (AB ° ~ 18.5 mag 
from the spectrum), comparable in brightness to the ob- 
servations presented here. The Wa ( ~ 15 A) reported in 
this paper and by Stickel et al. (|l993|) place B2 1308+326 
outside the traditional limits for BL Lacs, while previous 
observations place it well inside. 

The level of the continuum may be the most signifi- 
cant factor in changing the value of Wa for the Mg n line 
and not variation of the intrinsic line luminosity, since (a) 
both large Wa measurements were made near low luminos- 
ity states and (b) the line flux of 3.8 x 10 -15 ergcm~ 2 s _1 



reported by Miller et al. ( 1978|) is within a factor of two 



Table 4. Average radio, optical and x-ray luminosi- 
ties for samples of XBLs (Slew), RBLs (Uy) and FS- 
RQs. The Slew, 1 Jy and FSRQ sample data (Cols. 1, 
2 and 3 respectively) are from Fossati et al. ( [1998 ). Lu- 
minosities calculated from the contemporaneous data for 
B2 1308+326 are shown for comparison. The 5 GHz lu- 
minosity of B2 1308+326 is estimated from an extrapola- 
tion of the best-fit parabola to the contemporaneous data. 
The minimum recorded flux at 5 GHz and the maximum 
recorded flux in the V-band for B2 1308+326 are shown 
in parentheses. Luminosities are given in units of erg s . 



Band 


XBLs (Slew) 


RBLs (Uy) 


FSRQ 


B2 1308+326 


5 GHz 

V-band 
IkeV 


41.71 
44.91 
44.94 


43.69 
45.68 
44.72 


44.81 
46.58 
45.98 


45.4 (44.0) 

45.1 (46.7) 

44.9 



The existence of sources intermediate between BL Lac 
and QSOs is expected in the blazar unification paradigm 



presented by Fossati et al. (1998), where synchrotron peak 



frequency is correlated with source luminosity. X-ray, op- 
tical and radio luminosities for B2 1308+326 are compared 
to averages for the different blazar samples examined by 
Fossati et al. (1998) in Table 0. It is interesting to note 
that the luminosities for B2 1308+326 span the range be- 
tween the RBL-FSRQ averages. 

Sources intermediate between BL Lac and QSOs with 
luminous emission lines that arc SSC-dominatcd and have 



high synchrotron peak frequencies ( Perlman et al., 1998 
[Sambruna ct al., 2000 ) still pose serious proble ms to the 



unification scheme proposed by Fossati et al. (1998). It 
has recently been pointed out (Georganopoulos, 2000) 
that the observed differences between FSRQs with high 
synchrotron peak frequencies ( "blue quasars" ) and HBLs 
can be explained by taking into account the orientation 
of the jet as well as the intrinsic source luminosity. In 
a similar way earlier unification schemes have attempted 
to unify intrinsically lower luminosity FRI radio galax- 
ies with BL Lacs and the higher luminosity FRII radio 



of the value of 5.8 x 10~ 15 crgcm~ 2 s -1 of Stickel et al. galax ies with quasars ( [Padovani, 1992| phiscllini et al 



( |1993| ) |evcn though Wa differs substantially between the 
observations. No error was given for these measurements. 
Padovani (1992) has shown that in general the low values 
of Wa observed in BL Lacs are due to intrinsically less lu- 



1993). It now seems likely that quasars and BL Lacs form 



minous line fluxes than in quasars. Sambruna et al. (2000) 
find Mgn luminosities of Lj^gii ~(4-40)xl0 43 ergs _1 for 
a sample of three FSRQs, comparable to but higher than 
B2 1308+326, with L Mg „ = 1.5 x lO^ergs" 1 . It there- 
fore appears that the small values of Wa observed in 
B2 1308+326 are due mostly to the extremely bright con- 
tinuum. Unfortunately, it was not possible to obtain a 
line flux from these observations. Similar intermediate 
BL Lac-QSO behaviour is not unique to B2 1308+326, 
and has also been o bserved in PKS 0521-365 for example 
([Scarpa et al., 1999|) . 



two distinct populations of higher and lower luminosity 
AGN respectively, where the orientation of the jet influ- 
ences the luminosity and peak synchrotron frequency. The 
unusual intermediate properties of B2 1308+326 make it 
an ideal candidate for studying the overlap between the 
BL Lac and QSO categories. It may therefore represent 
the prototype of this kind of composite source. 

5.1. Gravitational microlensing 



The gravitational microlensing scenario (Gabuzda et al 



1993; McBrcen & Metcalfe, 1987), provides an alternate 



explanation for the BL Lac properties of B2 1308+326. Its 
high bolometric luminosity is comparable to that of other 
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gravitationally lensed candidates AO 0235+164 (Rabbette 
et al., 1 1996|) and PK S 0537-441 ( [Romero et al., 1999|; 



Pcrlman ct al., 1996]) and may be due to amplification 



by a foreground galaxy. Mgn absorption lines with very 
small W A (~ 0.4 A) were reported ( [Miller ct al., 1978| ) at 
z = 0.879 in B2 1308+326. These lines have not been con- 
firmed by later observations. Recent imaging with HST 



(Urry et al., 1999) did not confirm the detection of spa- 
tially extended emission ( Stickel et al., 1993| ) and no ev- 
idence was found for a foreground lensing galaxy. Such 
a galaxy is therefore limited to an I-band magnitude 
greater than 20.1, which, if the galaxy were at a red- 
shift z — 0.879 implies Mi > —23, thus not constrain- 
ing its luminosity significantly. The absorbing column to- 
wards B2 1308+326 is greater than that expected from 
our Galaxy (Fig. ||). The excess absorption could be due 
to an absorber at z < 0.997 or equally, to absorption in 
the host galaxy of B2 1308+326. The detection of excess 
absorption therefore does not exclude the possibility of a 
foreground lensing galaxy. A more detailed spectroscopic 
observation with XMM could determine the redshift of 
the foreground absorber. Sensitive, high-resolution optical 
spectroscopy could confirm the detection of Mg II absorp- 
tion lines (Miller et al., 1978) and thus also determine the 



redshift of such an absorber. 

5.2. Constraints on beaming in B2 1308+326 

Contemporaneous multiwavelength observations allow the 
determination of a lower limit on the doppler boost factor 
(8), the corresponding observed brightness temperature 
Tb, the maximum angle (</>) between the line of sight and 
the velocity of the relativistic jet as well as the minimum 
Lorentz factor of the jet, T. Ghisellini et al. ( 1995 ) have 
derived these values for B2 1308+326 with data from dif- 
ferent epochs, using the equations: 



5=(l + z)f(a)F n 



In {y b /v m ) 



F x e° d 



6 + 4q ,,a ,,5+3a 



1.77 x 10^(1 



v x v -iri 



F, 



J d u m 



K 



arctan 



2/3 a 



01 + 8 2 -l 



and 



28 



(1) 



(2) 



(3) 



(4) 



where F m and F x are the fluxes in Jy at the synchrotron 
self-absorption frequency, v m (GHz), and at the x-ray or 
optical energy, v x (keV) respectively; v b is the synchrotron 
high energy cut-off frequency; a is the spectral index of 
the thin synchrotron emission defined from F v oc v~ a ; a 
is the fastest reported apparent superluminal speed; 9^ is 



the VLBI FWHM core size flMarscher, 1977| ) in milliarc- 
seconds (mas) at v m \ z is the redshift of the source. A 
comparison of these values obtained at different epochs 
for b2 1308+326 is made in Table |. 

It has been noted that the largest source of error in 
Eq. p) is the deter mination of the VLBI core size, 6d 
(Madau et al., 1987). This value is crucial in determining 



the limit on S. 

Ghisellini et al. ( 1993| ) calculate 8 — 6.8 using 
0d = 0.5 mas at 5 GHz. High resolution VLBI images of 
B2 1308+326 were presented by Gabuzda et al. ( |1993| ) ob- 
tained in 1987 and 1989. Using the 5 GHz fluxes and cor- 
responding FWHM core sizes of Gabuzda et al. (1993), 8 
values of 122 and 9.4 were determined with the data pre- 
sented here. Clearly the core size is variable and further 
constraints on 8 would require VLBI imaging simultane- 
ous with observations at other wavelengths. 



Terasranta & Valtaoja (1994) use the radio variabil- 
ity timescales of B2 1308+326 to deduce a lower limit for 
the doppler boost factor of 8 = 1.9 and the correspond- 
ing maximum observed brightness temperature Te.max = 
6.4 x 10 12 K. The 22 GHz light curve used for this calcu- 
lation is plotted in Fig. 5. This method provides a useful 
independent measure of 8. A comparison of these values 
with those derived using Eqs. (p|) — (|4|) are given in Table 0. 
Two sets of values have been derived with the data ob- 
tained contemporaneously with ASCA u sing t he different 
core sizes reported by Gabuzda et al. ( 1993J ). The con- 
temporaneous data indicates that the synchrotron peak 
frequency is roughly 10 12 Hz (Fig. g) and this value was 
used for v\>. A similar fit to the highest archival fluxes in 
each waveband and the ROSAT data below the spectral 
break yields a peak synchrotron frequency of ~ 10 14 Hz. 
The V-band flux and frequency were used instead of the 
IkeV flux in these calculations, as the IkeV flux at the 
time of these observations was almost certainly dominated 
by inverse Compton emission. 

Three of the limits on 8 agree within a factor of five 
(Table g|). The extreme value of 5 = 122, obtained us- 
ing the core size 8d = 0.05 mas, is unlikely given the low 
flux state of the source during these observations. While 
the upper limits on (f> are consistent with expected jet 
angles for RBLs, they are also consistent with those of 
highly-polarised quasars (HPQs; Ghisellini et al. 1993, 
Urry & Padovani 1995). Quasars generally appear to have 
larger values of 8 than BL Lacs (Gabuzda et al., 1993; 
[Lcrasranta fc Valtaoja, 1994). A mean <5 of 8.3 for HPQs 



and 1.3 for BL Lacs (Gabuzda et al., 1993) indicates that 
B2 1308+326 with 8 = 9.4, is rather extreme for a typical 
BL Lac. 

6. Conclusions 

The blazar B2 1308+326 was observed contemporaneously 
at x-ray, optical and radio wavelengths in June 1996. The 
source was in a low synchrotron flux state at the time of 
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Table 5. Comparison of relativistic jet parameters calculated using contemporaneous data compared to those of 
Ghiscllini ct al. (1993) and Terasranta & Valtaoja (1994). F m is the flux at the synchrotron self-absorption frequency 
v mi F x is the x-ray or optical flux at the energy u x , fj, is the synchrotron high energy cut-off frequency, a is the spectral 
index of the thin synchrotron emission, (3 a is the highest superluminal speed, 0d is the VLBI FWHM core size at u m , 
z is the source redshift, S is a lower limit on the Doppler boost factor, Tb is the observed brightness temperature, <fi is 
an upper limit on the angle between the jet velocity and the line of sight and T is a lower limit on the Lorentz factor 
of the jet. 8, Tb, <f> and F were calculated using the data listed in that row. The third and fourth rows use data from 
observations made contemporaneously with the ASCA observations using the two different VLBI core sizes reported 



by Gabuzda et al. (1993) 



(Jy) (GHz) (Jy) 



(keV) 



(GHz) 



K' 



Vd 

(mas) 



T B <f> 

(10 12 K) (°) 



1.97 



1.22 c 
0.76 c 



Gh 



5.0 



Gh 



0.3 x 10 



-6 Gh 



5.0 
5.0 



Go 



Go 



0.18 x 10" 
0.18 x 10" 



,Gh 



2.28 x 10" 



1 x 10 5 Gh 0.75 G ' 1 



1 x W 0.80 

2.28 x 10~ 3 1 x 10 3 0.80 



15.75 
15.75 
15.75 
15.75 



0.5 



Gh 



0.05 
0.18 



Ga 



Ga 



0.996 

0.992 T 

0.997 

0.997 



Gh 



6.8 

1.9 T 

122 

9.4 



Gh 



LIT 
6.40 -7 
69.0 
3.35 



6.1 21.7 

7.2 68.0 
0.12 62.0 
5.37 18.0 



Data from Ghisellini et al. (1993) Ga Data from Gabuzda et al. (1995) 1 Data from Terasranta & Valtaoja (1994) 



observation and no variability was detected. The ROSAT 
data reveal an x-ray spectrum that is best fit by a broken 
power law with absorber model, with the break at 1.1 keV 
in the rest-frame of the blazar. The break is probably due 
to the emerging dominance of the IC over the synchrotron 
component. The frequency of the peak synchrotron emis- 
sion component appears to have decreased with decreasing 
flux, thus indicating that SSC may be the dominant emis- 
sion mechanism in this source. The x-ray IC component 
is unaffected by the large change in the synchrotron emis- 
sion. Mgn emission was observed with rest Wa ~15A, 
significantly different from the Wa values previously re- 
ported. The variable Wa emission in B2 1308+326 is prob- 
ably mostly due to the highly variable continuum. Al- 
though values of Wa are quite small, the line luminosity is 
only slightly lower than in a sample of QHBs. A lower limit 
on the Doppler boost factor obtained from the contempo- 
raneous data is consistent with expectations for HPQs but 
is higher than expected for BL Lacs. X-ray absorption at 
a level in excess of the Galactic value was detected and 
indicates the possible presence of a foreground absorber. 



B2 1308+326 is a typical RBL in terms of peak syn- 
chrotron power and optical variability, but with a very 
large bolometric luminosity, variable line emission and a 
high Doppler boost factor also has quasar-like attributes. 
It may be an intermediate or transitional AGN or a grav- 
itationally microlensed quasar, given its absorption above 
the Galactic level. Future x-ray observations could deter- 
mine the nature of the absorber. Further simultaneous 
multiwavelength observations of this possibly prototypi- 
cal intermediate blazar could determine the cause of both 
its BL Lac and its quasar-like properties, providing insight 
into the relation between BL Lacs and FSRQs. 



Acknowledgement. We would like to thank A. Celotti, 
for helpful comments that resulted in significant improve- 
ments to the paper. 



References 

Angel J. R. P., Stockman H. S., 1980, ARA&A 18, 321 
BevingtonP. R., Robinson D. K., 1992, In: Data reduction 

and error analysis for the physical sciences, 2nd ed., 

McGraw-Hill, New York 
Bottcher M., 1999, ApJ 515, L21 
Brown L. M. J., Robson E. I., Gear W. K. et al., 1989, 

ApJ 340, 129 
Catanese M., Bradbury S. M., Breslin A. C. et al., 1997, 

ApJ 487, L143 
Chiaberge M., Ghisellini G., 1999, MNRAS 306, 551 
Comastri A., Molendi S., Ghisellini G., 1995, MNRAS 

277, 297 
Costamante L., Ghisellini G., Giommi P. et al., 2000, 

astro-ph/0001410 
Dickey J. M., Lockman F. J., 1990, ARA&A 28, 215 
FanaroffB. L., Riley J. M., 1974, MNRAS 167, 31p 
Fiorucci M., Tosti G., Rizzi N., 1998, PASP 110, 105 
Fossati G., Celotti A., Ghisellini G., Maraschi L., 1997, 

MNRAS 289, 136 
Fossati G., Maraschi L., Celotti A., Comastri A., Ghisellini 

G., 1998, MNRAS 299, 433 
Fossati G., Celotti A., Chiaberge M. et al., 2000, astro- 

ph/0005066 
Gabuzda D. C, Kollgaard R. I., Roberts D. H., Wardle J. 

F. C, 1993, ApJ 410, 39 
Georganopoulos M., 2000, astro-ph/0008482 
Georganopoulos M., Marscher A. P., 1998a, ApJ 506, Lll 
Georganopoulos M., Marscher A. P., 1998b, ApJ 506, 621 
Ghisellini G., Maraschi L., 1989, ApJ 340, 181 



D. Watson et al.: ASCA and other contemporaneous observations of the blazar B2 1308+326 11 

Ghisellini G., Padovani P., Celotti A., Maraschi L., 1993, Rabbette M., McBreen B., Steel S., Smith N., 1996, A&A 

ApJ 407, 65 310, 1 

Ghisellini G., Celotti A., Fossati G., Maraschi L., Comas- Rabbette M., McBreen B., Smith N., Steel S., 1998, A&AS 

tri A., 1998, MNRAS 301, 451 129, 445 

Giommi P., Ansari S. G., Micol A., 1995, A&AS 109, 267 Romero G. E., Ccllonc S. A., Combi J. A., 1999, A&AS 
Giommi P., Massaro E., Chiappetti L. et al., 1999, A&A 135, 477 

351, 59 Sambruna R. M., 1997, ApJ 487, 536 

Gopal-Krishna, Subramanian K., 1991, Nat 349, 766 Sambruna R. M., Maraschi L., Urry C. M., 1996, ApJ 

Heidt J., Wagner S. J., 1998, A&A 329, 853 463, 444 

Inoue H., 1993, Experimental Astronomy (ISSN 0922- Sambruna R. M., Ghisellini G., Hooper E. et al., 1999, 

6435) 4, 1 ApJ 515, 140 

Iwasawa K., Fabian A. C, Nandra K., 1999, MNRAS 307, Sambruna R. M., Chou L. L., Urry C. M., 2000, ApJ 533, 

611 650 

Kollgaard R. I., Wardle J. F. C, Roberts D. H., 1990, AJ Scarpa R., Falomo R., 1997, A&A 325, 109 

100, 1057 Scarpa R., Urry C. M., Falomo R., Treves A., 1999, ApJ 

Kollgaard R. I., Wardle J. F. C, Roberts D. H., Gabuzda 526, 643 

D. C, 1992, AJ 104, 1687 Stickcl M., Fried J. W., Kiihr H., Padovani P., Urry C. M., 

Kollgaard R. I., Gabuzda D. C, Feigelson E. D., 1996, 1991, ApJ 374, 431 

ApJ 460, 174 Stickel M., Fried J. W., Kiihr H., 1993, A&AS 98, 393 

Kubo H., Takahashi T., Madejski G. et al., 1998, ApJ 504, Tanaka Y., Inoue H., Holt S. S., 1994, PASJ 46, L37 

693 Terasranta H., Valtaoja E., 1994, A&A 283, 51 

Lamer G., Brunncr H., Staubert R., 1996, A&A 311, 384 Terasranta H., Tornikoski M., Mujunen A. et al., 1998, 
Laurent-Muehleisen S. A., Kollgaard R. I., Feigelson E. D., A&AS 132, 305 

Brinkmann W., Siebert J., 1999, ApJ 525, 127 Urry C. M., Padovani P., 1995, PASP 107, 803 

Macomb D. J., Akcrlof C. W., Allcr H. D. ct al., 1995, Urry C. M., Falomo R., Scarpa R. ct al., 1999, ApJ 512, 

ApJ 449, L99 88 

Madau P., Ghisellini G., Persic M., 1987, MNRAS 224, Urry C. M., Scarpa R., O'Dowd M. et al., 2000, ApJ 532, 

257 816 

Maraschi L., Fossati G., Tavccchio F. et al., 1999, ApJ Watson D., Hanlon L., McBreen B. et al., 1999, A&A 345, 

526, L81 414 

Marcha M. J. M., Browne I. W. A., Impey C. D., Smith Wills B. J., Wills D., 1979, ApJ 41, 689 

P. S., 1996, MNRAS 281, 425 
Marscher A. P., 1977, ApJ 216, 244 
McBreen B., Metcalfe L., 1987, Nat 330, 348 
Miller J. S., French H. B., Hawley S. A., 1978, In: A. M. 

Wolfe (cd.), Pittsburgh Conference on BL Lac Objects, 

pp 176-191, University of Pittsburgh, Pittsburgh 
Morrison R., McCammon D., 1983, ApJ 270, 119 
Mufson S. L., Stein W. A., Wisniewski W. Z. et al., 1985, 

ApJ 288, 718 
Mukherjee R., Bottcher M., Hartman R. C. et al., 1999, 

ApJ 527, 132 
Nottale L., 1986, A&A 157, 383 
Ostriker J. P., Vietri M., 1990, Nat 344, 45 
Padovani P., 1992, MNRAS 257, 404 
Padovani P., Giommi P., 1995, ApJ 444, 567 
Pcrlman E. S., Stocke J. T., Schachter J. F. et al., 1996, 

ApJS 104, 251 
Perlman E. S., Padovani P., Giommi P. et al., 1998, AJ 

115, 1253 
Pfeffermann E., Briel U. G., 1986, Proceedings of SPIE 

597, 208 
Pian E., Falomo R., Ghisellini G. et al., 1996, ApJ 459, 

169 
Pian E., Vacanti G., Tagliaferri G. ct al., 1998, ApJ 492, 

L17 



